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ABSTRACT 

Spitzer IRAC observations of 15 metal-polluted white dwarfs reveal infrared excesses in the spectral 
energy distributions of HE 0110-5630, GD61, and HE 1349-2305. All three of these stars have 
helium-dominated atmospheres, and their infrared emissions are consistent with warm dust produced 
by the tidal destruction of (minor) planetary bodies. This study brings the number of metal-polluted, 
helium and hydrogen atmosphere white dwarfs surveyed with IRAC to 53 and 38 respectively. It also 
nearly doubles the number of metal-polluted helium-rich white dwarfs found to have closely orbiting 
dust by Spitzer. From the increased statistics for both atmospheric types with circumstellar dust, 
we derive a typical disk lifetime of log[tdisk(yr)] = 5.6 ± 1.1 (ranging from 3 x 10 4 — 5 x 10 6 yr). 
This assumes a relatively constant rate of accretion over the timescale where dust persists, which 
is uncertain. We find that the fraction of highly metal-polluted helium-rich white dwarfs that have 
an infrared excess detected by Spitzer is only 23 per cent, compared to 48 per cent for metal-polluted 
hydrogen-rich white dwarfs, and we conclude from this difference that the typical lifetime of dusty disks 
is somewhat shorter than the diffusion time scales of helium- rich white dwarf. We also find evidence for 
higher time-averaged accretion rates onto helium-rich stars compared to the instantaneous accretion 
rates onto hydrogen-rich stars; this is an indication that our picture of evolved star-planetary system 
interactions is incomplete. We discuss some speculative scenarios that can explain the observations. 
Subject headings: circumstellar matter — minor planets, asteroids: general — planetary systems — stars: 
abundances — white dwarfs 



1. INTRODUCTION 

Over the past decade it has become increasingly clear 
that planetary systems survive, at least in part, the late 
evolution of their host stars. The observational evi- 
dence supporting this conclusion comes from metal pollu- 
tion observed in white dwarf atmospheres, and the com- 
monly detecte d circumstellar disks o f solid a nd gaseous 
debris (e.s.. IZuckerman fc Becklinl Il987f 
IGansicke et all 120061. I2007L 12006 



been unambiguously identified as circumstellar dust 



Kilic et al 



Jura et al 



Ivon Hippel et all 120071 iBrinkworth et al.) 1 2009: 



Farihi et al.l2010bl : iDebes et al.ll2011t iFarihi et al.ll2011at 
Melis et alJl2011l Brinkworth et al. 2012 submitted, and 



a host of other references therein) . 

Due to the high surface gravities and limited radia- 
tive forces, heavy elements sink rapidly within the at- 
mospheres of cool (T c ff < 25, 000 K) white dwarf stars, 
resu lting in essent ially pure hydrogen or helium spec- 
tra (jKoesterl [20 09) . The downward diffusion of metals 
occurs over a timescale that is typically only days to 
years for hydrogen-rich (DA type) atmospheres (above 
~ 11,000K). Despite this, a notable fraction of cool 
DA white dwarfs exhibit abso rption features due to 
the presence of tra ce metals ([Zuckerman eiTaTI 120031 
IKoester et al.ll2005O ) , and these must be the r esult of ex- 
terna l sources and a sign of ongoing accretion ()Sion et al.1 
[1991 . 

In numerous cases, the origin of these pollutants have 
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Reach et al.1 120051 : 
20101^ Gaseous 



primarily via Spitzer studies (e.g. 
lJura et all [20091 IFarihi et all 1200a 
disks components are detected around a number of stars, 
and the analysis of their dynamics de monstrates that 
the material is l ocated within ~ 1R Q ( Gansick e et al.l 
I2006LI2007L 120081 ). where large solid bo dies should be de - 
stroyed by tidal gravitational forces (jDavidssonl I1999T) . 
The chemical abundances of these photospheric pollu- 
tants ca n be analytica lly linked to those of the accreted 
matter (|Koesterll2009f h In a number of cases, the closely 
orbiting dust has been confirmed to be silicate-rich, and 
typical of the material associated wit h planet formation 
(e.g. zodi acal a nd comet ary dust- iReach et al.1 120051 : 
Uura et al.l 120091: IReach et all 120091) . Studying plan- 
etary systems around white dwarfs therefore unlocks 
the potential to measure the bulk chemical composition 
of destroyed, and subsequently accreted, rocky plane- 
tary bodies such as asteroids, moons, or possibly major 



planets (e.g . IZuckerman et all [20071: iDufour et atl l2010l: 



planets (e .g. Zuckc-ma n et al.l ^00/: I 
IFarihi et afll2011al: Klein et al.ll201lh . 



The favored and successful model of a tida lly de- 
stroyed asteroid (|Graham et al.|[l990l : Uurall2003D is con- 



sistent with the observed disk properties (IGansicke et ahl 
20061: IFarihi et al.1 120101 iMelis et all 120111 : IDebes et al.l 
201 lh . the subsequent photospheric pollution, and the 
compositi on of both the orbiting and accreted mate- 
rial (e.g. ^^^^Ml^l^E^^^^^^l^l^). Such 
a catastrophic destruction is most readily achieved by 
a remnant planetary system with at least one major 
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planet ([Debes &: Sigurdssonl I2002D that perturbs a belt 
of smaller objects. Thus white dwarfs with disks and 
photospheric metals may harbor complex planetary sys- 
tems. 

A signific ant uncertainty is the typica l lifetime of the 
dust disks (| Jural 120081 iKilic et al.1 (20081 ) . yet this is an 
important indicator of the mass of the parent body (or 
bodies) that generated the observed debris. There are 
numerous stars with cooling ages greater than 1 Gyr 
that exhibit atmo spheric metal pollution, such as the 
prototype, y Ma 2 (Ivan Maanenl Il917t iGreensteml 119561 : 
lWeidemannlll960l ) . Koester et all 1)20111 ) recently identi- 
fied 26 white dwarfs with temperatures 5000 — 8000 K 
and 10 20 — 10 23 g of accreted metals in their convec- 
tion zones. These stars have typical (log g = 8) cooling 
ages of 1 — 6 Gyr but must have accreted material re- 
cently, within the past few million years. Interestingly, 
there is only a single (and anomalous) infrared excess 
aro und stars with co oling ages older than 1 Gyr (G166- 
58: lFarihi et~aTll2008fl ■ so this may be viewed as an upper 
limit for typical disk lifetimes. 

From theoretical considerations, Rafiko3 (|2011a| ) finds 
the lifetimes of the compact dust disks around white 
dwarfs should be of the order 10 6 yr, when domi- 
nated by Poynting-Robertson drag. However, Poynting- 
Robertson drag cannot produce the highest (average) ac- 
cretion rates inferred for the helium atmosphere stars 
(DB) with metals, which are on the order of 10 10 — 
^gs" 1 dFarihi et all I2010bl ) IXu fc Jural ([2011) ex- 
tend the estimates of iRafikovl ([201 lal ) and find a factor 
of five higher accretion rates, but this still does not lead 
to high enough rates. To produce higher accretion rates, 
gas resulting from subli mated du s t, produced at the inner 
edge of the dust disk (|Jurall2008l: iFarihi et al.ll2009n , effi- 
ciently transports angular moment um outw a rd and fuels 
a more rapid in-fall of material. IRafikovl (|2011bD cal- 
culates that the lifetime of a 10 22 g disk is reduced to 
several 10 4 yr. 

We report the results of a combined Spitzer IRAC 
Cycle 5 and 6 search for disks around 15 cool, metal- 
contaminated white dwarfs. Three of the observed 
stars show infrared excess that is the trademark of 
warm dust orbiting within the tidal disruption radius: 
HE 0110-5630, GD61, and HE 1349-2305. These three 
stars significantly increase the number of metal-polluted 
DB white dwarfs with disks and enable a statistical es- 
timate of a typical disk lifetime based on comparison 
with their DA counterparts. The observations are de- 
scribed in §[2J the analysis and disk modeling in §[3] In 
§2] the statistics of dusty white dwarfs is revisited, and 
our method for estimating disk lifetimes is described. 

2. OBSERVATIONS 

We were awarded Spitzer Space Telescope time in Cy- 
cle 5 (program 50340) and Cycle 6 (program 60119), 
to search for infrared excesses from circumstellar dust 
around 16 white dwarfs published as metal-enriched (Ta- 
ble[T]). These targets were chosen from the literature, and 
significantly increase the number of DBZ (for simplicity, 
we use 'DBZ' to refer to all helium-rich subtypes with 
metals, including those with trace hydrogen) stars ob- 
served with Spitzer IRAC, bringing their numbers onto 
par with those of the DAZ stars. One of the targets, 
PHL131, was subsequently found by one of us (DK) to 



be a 40, 000 K white dwarf with interstellar calcium ab- 
sorption, and thus we report its fluxes, but exclude it 
from the rest of the study. 

Imaging observations were obtained for each of the 
white d warf targets usin g the Infrared Array Camera 
(IRAC: Fazio et all 120041 ) . Exposures were taken in the 

4.5 and 7.9 /im channels using the medium-scale, cycling 
dither pattern for the Cycle 5 targets objects and in the 

3.6 and 4.5 /im channels for the Cycle 6 objects. For 
GD61 and NLTT 51844, images were also obtained with 
the blue Peak- Up Imaging array of the Infrared Spectro- 
graph at 15.6 /im, this time using the small-scale, cycling 
dither pattern. 

Three of the target whit e dwarfs were cho sen from 
the Hamburg Schmidt (HS; lHagen et allll995l ) and the 
Hamburg European S outhern Observatory Schmidt (HE; 
Wisotzki et"aLl I1996D quasar surveys. As discussed in 
Farihi et all (|2010b[ ). the SIMBAD coordinates for these 
objects are often inaccurate by up to a few arc minutes, 
and Table[2]gives correct positions for these sources based 
on their IRAC images. 

2.1. Data Analysis 

The IRAC image photometry was performed on the 
individual BCD frames downloaded from the Spitzer 
archive and reduced using SI 8. 7.0. These were cor- 
rected for array-location dependence, as described in 
the IRAC Data HandboolQ. We used the point source 
extraction package a pex MULTIFRAME within MOPEX 
(|Makovoz et al .1120061 ) to perform PSF fitting photome- 
try on the individual BCDs. We do not apply a color 
correction because we use the isophotal wavelengths in 
our analysis. The magnitude of the correction is there- 
fore negligible compared to our uncertainties. The pixel 
phase correction for channel 1 data is also minimal in 
comparison to the uncertainties and therefore was not 
applied. A minimum 5 per cent uncertainty on the flux 
density was assigned to ac count for the syst ematic un- 
certainties in this method ([Reach et al.ll2005|). 

The Infrared Spectrometer (IRS iHouck et all l2004f ) 
Peak-up Imaging was analyzed according to the method 
of Brinkworth et al. (2011, submitted). The sky back- 
ground was found to be variable over the array (see Fig- 
ure!]]), so the individual BCDs were initially median- 
combined to make both a flat and sky, using 3<r rejec- 
tion. This was then scaled to the overall median of the 
Post-BCD mosaic and subtracted from each of the in- 
dividual BCDs. To completely reduce the sky level to 
zero, the edges of the image were temporarily discarded 
and a median of the center of the image was found and 
removed. The sky can therefore assumed to be zero and 
an infinite sky annulus is effectively used during photom- 
etry. The flat-fielded and median-subtracted BCDs were 
mosaicked using MOPEX with a dual-outlier rejection. 
The pixel scale for the mosaic was set to the default l'/8 
pixel -1 . The photometry was performed with APPHOT 
within IRAF using a 3 pixel aperture radius, but no sky 
subtraction. The aperture correction to the calibration 
aperture sizes, as provided by the Spitzer Science Cen- 
ter, was performed and the fluxes were converted from 
MJy sr _1 to mJy. As an estimate of the error, a series 
of other apertures were placed around the target on the 

1 See http://ssc.spitzer.caltech.edu/irac/dh 
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mosaic and the standard deviation of fluxes was taken as 
the uncertainty. All measured Spitzer fluxes are listed in 
Tabled 

2.2. Near-Infrared Observations 

Supplemental near-infrared photometry for most tar- 
get stars was obtained on the 23 March 2011 with 
the William Herschel Telescope using the Long-Slit In- 
termediate Resolution Infrared Spectrograph (LIRIS; 
iManchado et al.lll998t ). and on 10-12 August 2011 with 
the New Technology T elescope using Son of Isaac (SOFI; 
iMoorwood et al.lfl998H . Images taken in a 9-point dither 
pattern were obtained the J, H, and i^ s -band filters with 
typical total exposure times of 270 s in clear conditions; 
3 standard star fields were observed in a similar man- 
ner for photometric zero-point calibration. These data 
were reduced in the standard manner, by subtracting a 
median sky from each image in the dithered stack, flat 
fielding (using sky flats), then averaging and recombining 
frames. 

LIRIS suffers from what is known as a detector re- 
set anomaly, which appears in certain frames as a dis- 
continuous jump (in dark current) between the upper 
and the lower two quadrants. To remove this unwanted 
signal, after flat fielding and sky subtraction, the de- 
tector rows were collapsed into a median column (with 
real sources rejected), and subsequently subtracted from 
the entire two dimensional image. The resulting fully- 
reduced frames exhibit smooth backgrounds, free of the 
anomalous gradient. Analogously, SOFI exhibits signifi- 
cant cross-talk between quadrants, which becomes appar- 
ent for bright sources. The IRAF script CROSSTALK was 
employed on each of the raw frames prior to processing 
as described above, effectively removing any unwanted 
artifacts. 

Aperture photometry of standard stars and relatively 
bright targets was performed using r 4" aperture radii 
and sky annuli of ranging between 5" and 8". For rel- 
atively faint targets or those with neighboring sources, 
smaller apertures were employed with corrections de- 
rived from several brighter stars within the same image 
field and filter. In a few cases, PSF-fitting photometry 
(i.e., daophot) was used as a second method in addition 
to photometry with small apertures. All data taken in 
the K s -band filter were flux-calibrated using ARNICA 
(|Hunt et al.l 119981 ) if-band standard star photometry. 
The measured photometry can be found in TablelU 

3. ANALYSIS AND RESULTS 

3.1. Removal of Nearby Background Source Flux 

The IRAC mosaicked images of GD61, CBS 127, and 
HE 1349—2305 all reveal nearby background objects that 
have the potential to contaminate aperture photometry 
(see Figured] and Figure[3] for the images of CD 61 and 
HE 1349— 2305 respectively). The background object 
near GD 61 is separated by l'/9 at P.A. 25° in the images 
taken with LIRIS (epoch 2011.2), and is readily resolved 
in the H and K bands. The near-infrared brightness 
of this neighbor is (J,H,K) = (17.7,17.0,16.7) mag and 
never approaches the apparent magnitude of the white 
dwarf. Based on this fact and the colors of the neighbor, 
it must be a background object. This source contributes 
weakly in the 4.5 /im image, where it is well separated 



by apex, and is no longer seen at 7.9 fim. CBS 127 and 
HE 1349-2305 both have neighbors 2'.'8 distant (north 
and west respectively); these are satisfactorily separated 
by PSF fitting with apex. 

3.2. Spectral Energy Distributions 

FiguresHHH illustrate the spectral energy distributions 
(SEDs) of the 15 metal-rich white dwarfs listed in Ta- 
bled In addition to the supplementary near-infrared 
photometry described above, we include shorter wave- 
length photometry from a va riety of literature a nd cat- 
alog s ources (e.g., GALEX: [Martin et all [20051 SPSS 



DR7; lAbazaiian et all 120091 DENIS; lEpchtein et al 
1998, and CMC; [Copenhagen University Obs. et al 
2006). In cases where independent JHK observations 



were not taken for th is program, dat a from 2MASS 
(ISkrutskie et al.l I2006D and UKIDSS (|Lawrence et al.l 
120071) were used. Where possible, a comparison of the 
IRAC photometry was made to data available at simi- 
lar wavelengt hs from the Wide-fi eld Infrared Survey Ex- 
plorer ( wise [MiiMiiEiniloia) . 

For flux scaling of the spectral models in the figures, 
the most reliable photometry was used, which was pri- 
marily SDSS ugriz and the JHK data taken specifically 
for this program. Specifically, the GALEX fluxes were 
not used to constrain the white dwarf atmospheric mod- 
els because 1) metal lines may significantly suppress these 
fluxes compared with a pure helium (or hydrogen) model 
atmosphere and 2) the interstellar reddening to these ob- 
jects is not well known. 

Tabled] lists the best available stellar parameters from 
the literature for each of the science targets. We used 
these parameters as an initial estimate for the white 
dwarf spectral models. A grid of s y nthet ic white dwarf 
atmospheric spectra from iKoesterl ([20101 ) was interpo- 
lated to these temperature and surface gravity estimates. 
After comparing the fluxes from the models with the 
most reliable optical and near-infrared data, in some 
cases we found that a slightly different effective temper- 
ature provided a superior fit to the data (shown in the 
upper right corner of the relevant figure). These param- 
eters thus used in the SED fits do not reflect an indepen- 
dent parameter determination, but rather our attempt 
to best constrain the infrared photospheric emission. 

3.3. Stars With An Infrared Excess 

We model stars with an infrared excess as a combi- 
nation of the white dwarf photosphere and an optically 
thick, geometrically thin disk with a temperature pro- 



file T disk oc r^, with = -3/4 ([Adams et al.l fl988l: 



IChiang fc Goldreichl[l997t [Jurl|2003D. The ratio of stel- 
lar radius to distance sets the absolute scale of the white 
dwarf photospheric flux, and is fixed in our modeling 
based on the best available parameters. The free param- 
eters in the disk model are the inner disk radius (R- m ), 
the outer disk radius (i? ut) and the disk inclination (i). 
A grid of disk models was calculated with inner (T; n ) 
and outer disk temperatures (T out ; each corresponding 
to a value of Ri n and i? ou t) ranging from 100 to 1800 K 
in steps of 50 K, and inclination ranging from to 90°in 
steps of 5°. T out was also fixed to be cooler than T m . A 
least x 2 method was used to fit the H-, K- and IRAC- 
band infrared fluxes and provide an estimate of the un- 
certainties. 
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In Figure[5j the panels show a slice through the x 2 cube 
at the best fit solution. The upper panel of each pair dis- 
plays a slice at the best fitting T out , whilst the best fit 
inclination defines the slice in the lower panel. Regions of 
high x 2 are shown as dark areas and the least \ 2 solution 
is marked as a red circle and its flux is displayed in Fig- 
ured! Solid black lines show the la contours around the 
minimum, while 2 and 3<r contours are shown as dashed 
gray lines. The upper left of each lower panel is excluded 
by the condition T; n > T out . 

The flat disk models naturally possess a modest de- 
gree of uncertainty. The three free parameters (2j n , T out 
and i) are somewhat degenerate in how they combine 
to determine the viewed solid angle of the disk. Even 
in the case when longer wavelength data (such as 24 /im 
phot ometry) are ava ilable, this degeneracy is not bro- 
ken ()Jura et al.H20Q7D . Broadly speaking, the inner disk 
edges are fairly well constrained by their 2.2 and 3.6 /im 
emission. The 4.5 to 7.9 /im flux, however, can be re- 
produced by a relatively wide temperature range and a 
higher inclination, or a narrower temperature range and a 
lower inclination. Importantly, the newly detected disks 
mimic the emission seen at more than one dozen dusty 
whit e dwarfs observed with Spitzer IRAC and MIPS 
(e.g. iFarihi et al.l I2009D whose outer disk radii are con- 
sistently within the Roche limit d > 1 km solid bodies 
(~ 1R©). This corresponds to an outer disk tempera- 
ture near 500 K, which is a likely lower limit for T out when 
interpreting the x 2 surfaces. All together, the infrared 
emissions are precisely that expected for dust particles 
resulting from the tidal destruction of a large planetary 
body, whose constituent elements now rain onto and pol- 
lute the surfaces of the host stars. 

HE 0110-5630. The available photometry for 
HE 0110-5630 consists of GALEX far- and near-UV, 
DENIS /-band and J, H and K measurements taken 
with SOFI. The photospheric emission is relatively well 
constrained by the J- and -ff-band fluxes and, therefore, 
so is the extrapolation to longer infrared wavelengths. 
The measured IRAC excess for HE 0110— 5630 is consis- 
tent with warm circumstellar dust and is shown in Fig- 
ured! The best fit disk model has an inner and outer 
temperature 1,000 and 900 K, respectively, at an incli- 
nation of i = 60°. This corresponds to a narrow ring, 
however, as can be seen in the x 2 panels, the solution is 
degenerate with respect to inner disk temperatures above 
800 K, and inclination, and the outer temperature is not 
at all constrained. The y 2 surfaces shown in FigureEl 
show that a vast range in disk temperatures could pro- 
vide a fit to the excess. We therefore cannot confidently 
estimate any of the disk parameters. 

GD61. From analysis of IRAC images and Peak- Up 
Imaging, th i s objec t was found to have an excess in 
IFarihi et al.l ()2011aD . Here we include the JHK data 
taken with LIRIS to better constrain the stellar flux. 
These images resolve a nearby (background) source (Fig- 
ure^ that almost certainly contaminates the 2MASS 
photometry used in previous fits. From its spectral en- 
ergy distribution and brightness, this object is not asso- 
ciated with the white dwarf. In the Spitzer IRAC im- 
ages, the PSF of the white dwarf and the contaminant 
overlap slightly, however, the apex PSF-fitting photom- 
etry cleanly separates the two objects. We can therefore 
confirm the presence of a disk from the IRAC and IRS 



fluxes. 

When performing \ 2 fitting of the infrared excess, we 
exclude the 15.6 flux due to a possible increase owing 
to silicate emission, and therefore the disk continuum 
emission would be overestimated. Such emission has 
been detected in all eight dusty white dwarfs observed 
spect roscopically with IRS (|Reach et al.ll2005t Uura et all 
2009). A disk model with inner and outer-disk tempera- 
tures of 1, 450 and 800 K, respectively, and an inclination 
i = 85° is the best fit to the excess (Figured] and Fig- 
ure[5]). The disk model temperature and inclination are 
again degenerate, and a large range of parameters would 
fit the infrared fluxes. We therefore find the values re- 
ported p reviously (r in = 1, 300K, T out = 1,000K and 
i = 79°; Fari hi et al.l 12011a) are not inconsistent with 
ours, given the uncertainties in temperature and inclina- 
tion. 

HE 1349-2305. After deconvolving the white dwarf 
and neighboring, background source f§ l3.1l) . a significant 
infrared excess is still found over the stellar model (Fig- 
ured]). The white dwarf photosphere is well constrained 
by short wavelength and near-infrared photometry shown 
in Table[H WISE 3.4 and 4.6 /im fluxes are also available 
for this object, however, the nearby source seen in Fig- 
ure^ is within the WISE beam width and therefore sig- 
nificantly contaminates the WISE fluxes. We therefore 
do not include these fluxes in the disk modeling. The 
excess is best fitted with a disk having T- m = 1,700K, 
Tout = 550 K and an inclination of i — 85°. The uncer- 
tainties are rather modest in this case. 

NLTT 51844- This was one of two white dwarfs bright 
enough to attempt blue Peak-Up image using IRS. The 
15.6 /im flux measurement is in excess over the white 
dwarf model by slightly more than 3cr, similar to the 
MIPS 24 /im photometric exces ses detected at G238-44 
and G180-57 (|Farihi et all [2009). Spitzer galaxy counts 
at these wavelengths are 6 x 10~ 4 galaxies per square arc 
second (jMarleau et al.l I2004T ) for the brightness needed 
to produce an apparent excess in these cases. Therefore 
the probability that such a source falls within a radius 
of two IRS Blue Peak-Up FWHM (3'.'8) can be as high 
as 10 per cent. Thus the probability that three sources 
have excesses at roughly 20 /im (and not at shorter wave- 
lengths) due to background objects, out of 30 stars sur- 
veyed at these longer wavelengths, can be as high as 
65 per cent. 

While we cannot yet rule out a chance alignment in 
these cases, the astrometric position of NLTT 51844 in 
the IRS image mosaic coincides within 2" of its position 
in the IRAC images, based on the absolute Spitzer as- 
trometry provided in the processed image files. Both 
NLTT 51844 and G180-57 are cool, helium-rich stars 
where metals can be resident for Myr timescales, and 
hence their accretion history is unconstrained. In any 
case, dust that emits only at such long wavelengths must 
be sufficiently far from the white dwarf that it cannot be 
the immediate source of the atmospheric metals. Still, if 
the data for NLTT 51844 (as well as G238-44 and G180- 
57) are not spurious, we may be seeing an important clue 
to the nature of numerous polluted white dwarfs without 
obvious infrared excesses. 
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3.4. Stars Without An Infrared Excess 

The remaining white dwarfs (Figures^ and [7]) all show 
no evidence for infrared excess emission. They are rel- 
atively well constrained in the optical and near-infrared 
and therefore we can extrapolate with some certainty to 
IRAC wavelengths. Below we provide notes on a few 
particular stars. 

HE 0446— 2531. This star is perhaps the most highly 
polluted white dwarf observed by Spitzer and yet 
does not have an in frared excess. If the analysis of 
iFriedrich et al.l (|2000D is correct, this star contains 8 x 
10 24 g of metals in its atmosphere - a mass approach- 
ing that of Pluto. The long diffusion timescale for 
heavy elements in this star allows for the possibility 
that a disk has been fully dissipated, but the very large 
mass involved suggests we are still within at most a 
few diffusion timescales of the accretion event. Two 
possibilities are: 1) accretion is ongoing from gaseous 
debris, or 2) accretion has ended but total mass in- 
volved was at least an order of magnitude higher. It 
is noteworthy th at this star and HS 22 53+8023 are both 
highly polluted (jFriedrich et al.l [20001 ) yet lack infrared 
excess, with HE 0446— 2531 having nearly ten times 
more meta l mass in its convection zone. However, 
iKlein et al.1 (|2011| ) found that HS 2253+8023 has an or- 
der of_jnagnitude less_calcium than originally reported 
by IFriedrich et all (|2000D . and therefore the same may 
be true of HE 0446-2531. 

HE 0449—2554- A partly-resolved extension is seen in 
the IRAC 7.9 /im image, and the combined source is off- 
set slightly from the position of the white dwarf in the 
4.5 /jm image. While the SOFI JHK images of this star 
do not reveal any sources within a few arc seconds of the 
science target, the marginal, apparent excess seen in the 
SED (Figure[6|) is probably caused by a background ob- 
ject. Better data are needed to rule out contamination 
in IRAC. 

CBS 127. As mentioned above, a background object 
lies within a few arc seconds of the white dwarf, and could 
potentially contaminate moderate-size apertures used for 
IRAC photometry. We used PSF-fitting routines within 
APEX to ensure both stars were photometrically disen- 
tangled. There is no infrared excess at this white dwarf. 

HE 1350— 1612 was observed in both the Cycle 5 and 6 
programs and images were thus taken at 3.6, 4.5 (twice), 
and 7.9 /xm. The average 4.5 /xm flux is shown in TableQ] 
and Figure[7l The white dwarf was not detected in the 
7.9 /im image, and thus a 3a upper limit is given. Nev- 
ertheless, all the infrared data are consistent with the 
predicted white dwarf photosphere. 

PHL131 is a 40,000 K white dwarf with interstellar 
calcium absorption, and hence not metal-rich. We do 
not analyze this star but its IRAC fluxes are reported in 
Tabled 

4. DISCUSSION 

4.1. Updated Statistics and Accretion Rates 

Of the 15 metal-polluted stars surveyed, three have an 
infrared excess consistent with circumstellar dust disk, 
and indicative of a remnant planetary system. Includ- 
ing the objects presented here, the total number of 
metal-polluted white dwarfs observed with Spitzer in 
Cycle 1 through 7 is 92. This includes the 52 objects 



from Table 4 of iFarihi et al.l (|2009l ). 8 from Table 1 of 



Farihi et al.l (l2010bl). 3 white dwa rfs with gas disks from 
Gansicke et all (120061 [20071 [2008] ) . 14 DBZ white dwarfs 



from lXu fe Jural (|2012f ) and the 15 objects studied here. 
Of these 92 objects, DAZ and DBZ-type white dwarfs 
are represented in proportions of 39 and 53 respectively. 

O ur expanded study corroborates the prev ious find- 
ing (|Kilic et al.1120081: IFarihi et al.ll2009l l2010bD that the 
detection of infrared excess is less frequent among DBZ 
white dwarfs than among their DAZ counterparts. 11 
of 38 (29 per cent) surveyed DAZ stars have circumstel- 
lar dust, whereas the same fraction for surveyed DBZ- 
type white dwarfs is only 7 of 52 (13 per cent j3- Be- 
cause one would not expect disks to preferentially form 
around stars based on their atmospheric properties, the 
difference must be due to the longer diffusion timescales 
(|Paauette et all [1981 iKoester fcWilkenl 120061) . We will 
return to this point in more detail below. 

Figurc[5] illustrates the time-averaged metal accretion 
rate versus cooling age for the 88 metal-rich white 
dwarfs now obse rved by Spitzer with published abun- 
dances (four of IXu fc Jural 120121 stars, including one 
with a detected infrared excess, do not ha ve published 
abunda nces) , updating Figure 10 from IFarihi et al.l 
(|2010bf) . All accretion rates were calculated based on 
the accumulated metal a bundances, using Equation 2 of 
IKoester fc Wilkenl (|2006l ). as performed in IFarihi et al.l 
(I2009D. One diffe rence between the method used in 
Farihi et al.l ((2009) and here is that in previous stud- 
ies instantaneous and time-averaged accretion rates were 
calculated assuming the in-falli ng material either had So - 
lar calcium abundances (1 : 43: IKoester fc Wil kcn 2006), 
or 1 p er cent Solar (1 : 109; i.e. metals only; Uura et all 
l2007t IFarihi eFall [2001) by mass fraction, based on the 
measured calcium abundance. 

While the latter approach is still likely to be broadly 
correct, recent progress in the field enables a more ac- 
curate estimation. There are now of order 10 metal- 
polluted whi te dwarfs with measured Mg, Si, Ca. and F e 
abundances (jZuckerman et alj 120031 . 120071 I2010L 120111 
and a handful also with O abundances (Dufour et al.l 



2010; 



2011 



Klein et al. 


2010 


Mclis et al. 


2011 



Vennes et al. 


2010; 


Venncs et al. 


2011) 



Based on 



these data, it is increasingly clear that the accreted ma- 
terial has a composition similar to t hat of rocky, ter- 
restria l material of the Solar System. IZuckerman et akl 
(2010) find that Ca represents, on average, close to 1 
part in 60, by mass, of all the accreted heavy elements. 
For comparison, thi s ratio is 1 part in 62.5 for the Earth 
( Alleg re et al.lll995T ). Based on these facts, Figure[5] and 
Table[T] employ accretion rates, and convective envelope 
masses, assuming Ca is 1.6percent (1/62.5) of the total 
mass. To convert between the accretion rates quoted in 
IFarihi et all (|2009l I2010bl) and those shown in Figure^ 
one must multiply by a factor of 109/62.5 for DAZ white 
dwarfs and 43/62.5 for DBZ white dwarfs. 

Compared to the same Figures shown in previous (re- 
lated) papers, the accretion rates for the DAZ white 

2 We exclude here the two DAZ and one DBZ tha t were already 
known to have gaseous disks (Gansicke ct al. 2008) p rior to their 



Spitzer observations. The stars discussed in Farihi ct al. (2011b) 
were also not included because they were observed with Spitzer 
after they were found to have an infrared excess in lUirven et al.l 
l20Tll ). 



6 



Girven et al. 



dwarfs are therefore slightly decreased, while their DBZ 
counterparts are slightly increased. 

4.2. A Simple Estimate of the Disk Lifetime 

The mass of metals in the convective envelope of DBZ- 
type white dwarfs provides a lower limit on the total mass 
of the destroyed parent body. From the DAZ stars, we 
can infer the instantaneous accretion rates because one 
can safely assume accretion-diffusion equilibrium. There 
is no a priori reason that disks aro und hydrogen an d 
helium-rich stars should be different ([Jura et~aLl[2007l) . 
Therefore, assuming the inferred accretion rates of DAZ 
stars represent an accurate cross section over disk life- 
times, and similarly for the observed metal masses con- 
tained within DBZ stellar envelopes, we can obtain a 
typical disk lifetime by combining the average of these 
two observed quantities, i.e. 

I ^ (Mz) DBZ, disk 

(dM z /dt) DAZ,disk 

Figure|9] demonstrates that the mass of metals in the 
convective envelope of DBZ white dwarfs is not a strong 
function of temperature, and hence cooling age. There- 
fore, calculating the average convective envelope metal 
mass across cooling ages should be robust. In contrast, 
the convective envelope depth increas es substantiall y in 
this temperature range for DAZ stars (jKocstcr 200jJ , as 
can be seen by the large increase in convective envelope 
metal masses at longer cooling ages in FigureO 

Among 13 DAZ white dwarfs observed to have infrared 
excess with Spitzer, the average metal accretion rate is 
9.7 x 10 8 gs _1 . Among 8 DBZ-type stars with dust de- 
tected by Spitzei^ the average metal content of the con- 
vection zone is 4.1 x 10 22 g, and we thus estimate idisk — 
1.3 x 10 6 yr. However, Figured] shows that the accretion 
rates and convection zone metal content vary by orders of 
magnitude, and suggests that using the logarithmic mean 
and standard deviation, (log[dM z /dt (gs -1 )]} = 8.8±0.4 
and (log[M z (g)]) = 21.9 ± 1.1, may be a more appropri- 
ate choice. The resulting estimate for the disk lifetime is 
log[idisk (yr)] = 5.6 ± 1.1. Because the metal masses in 
the DBZ stars are lower limits on the total parent body 
masses, it is reasonable to presume the same for the disk 
lifetimes. 

At least 1 per cent of all white dwar fs with cooling ages 
less than 0.5 Gyr have dusty disks (iFarihi et al.l 120091 : 
I Girven et all 120111 : iSteele et afl 1201 lh . Assuming that 
all white dwarfs host (remnants of) planetary systems 
and go through intermittent phases of debris, one would 
expect any given star to exhibit detectable infrared ex- 
cess for 0.01 x 0.5 x 10 9 yr = 5 x 10 6 yr, which is broadly 
consistent with our estimate above. 

4.3. How Complex is the Evolution of the Dust Disks? 

Different possible scenarios for the evolution of dust 
disks, and corresponding estimates for t heir l i fetime s 
have been discussed in more detail bv [Jural (l200ll . 
IRafikovi HonaE), and IBochkarev fc Rafikovl ([201 lh . 



A main motivation of the work of Uural (|2008f ) was 
to explain the existence of metal-polluted white dwarfs 
without infrared excess detection by the continuous ac- 
cretion from a gaseous disk that is replenished by the 
repeated tidal destruction of multiple small asteroids. In 
this scenario, dust disks are associated w ith the disrup- 
tion of massive asteroids, and I Jural (|2008f) estimated the 
lifetimes of the gaseous and dust disks to be ~ 0.5 x 10 4 yr 
and ~ 1.5 x 10 5 yr, respectively, though he underlined 
the uncertainty of these estimates. I Jural ((2008) discussed 
his scenario in the context of a much smaller sample of 
metal-polluted white dwarfs, in which the majority of 
high dM/dt stars were displaying infrared excess iden- 
tifying the presence of dusty disks. It is unlikely that 
the much larger accretion rates that are now known for 
some stars without infrared excess can be explained by 
multiple impacts of small asteroids. 

Rafikov (2011a) showed that Poynting- Robertson drag 
can explain accretion rates of ~ 10 8 g/s, and estimated 
the lifetime of massive (~ 10 22 g) dust disks to several 
Myr. The higher accretion rates observed in a num- 
ber of white dwarfs requir e the pre s ence of gas that in- 
creases the viscosity, and Rafikov (2011b) shows that 
sublimation of the inner disk may be sufficient to lead 
to a runaway evolution which produces accretion rates 
of 10 10 — 10 11 g/s, exhausting the disk in ~ 10 5 yr. 
IRafikovi (|2011bf ) explicitly excluded the generation of gas 
by multiple impacts of additional asteroids. Extending 
the analysis to the global e volution of the debris disks, 
IBochkarev fc Rafikovl (|201lD confirm that massive disks 
have expected life times of several Myr if subject to 
Poynting- Robertson drag only, but that sublimation may 
substantially shorten these estimates. 

Our disk lifetime estimate in the previous section is 
broadly compatible with all the estimates summarized 
above. However, our approach assumed that accretion 
rates do not change significantly as a function of disk 
age, and that the Spitzer sample is rep resentative of the 
poten tial changes in the accretion rate. [Zuckerma n et al.1 
( 2010) suggested that accretion onto a white dwarf could 
proceed through three distinct phases; a build-up, a 
steady state, and a decline, where the steady state life- 
time should be orders of magnitude longer than the be- 
ginning and ending episodes. If correct, and disk life- 
times are typically ~ 10 6 yr, then it will be very unlikely 
for an arbitrary white dwarf with circumstellar dust to 
be observed in anything other than a steady state phase. 

Inspection of Figure[8] shows, as already mentioned in 
§ 14. 1[ that the fraction of DBZ with Spitzer infrared 
excess is substantially smaller than that of DAZ. The 
detection of dust seems to be associated with a mini- 
mum accretion rate, and if we only consider white dwarfs 
with dM/dt > 10 8 gs _1 , we find that 10/21 (48 per cent) 
of DAZ have infrared excess, compared to only 7/30 
(23 per cent) for the DBzQ This difference is, however, 
subject to small number statistics and future observa- 
tions may show that both DAZ and DBZ are equally 
likely to show an infrared excess. If there is a real dis- 
parity, this fact is most likely related to the one signif- 
icant difference between DBZ and DAZ: their diffusion 



3 This does not in c lude SPSS J220934.84+122336.5, the DBZ 
star from IXu &~J ura ( 2012) with an infrared excess, because no 
calcium abundances are published. 



4 Again excluding the two D AZ and one DBZ that were identified 
because of their gaseous discs (Gansickc et al. 2008) 
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time scales, and allows some independent gauge of the 
disk lifetime. Assuming a typical lifetime of the disks 
that is similar to or greater than the diffusion time scale 
of the DBZ stars {t disk > t diS ~ 10 5 - 10 6 ), we would 
expect a similar ratio of DAZ and DBZ stars exhibiting 
infrared excess. By contrast, if idisk "C tdiff , it would be 
very unlikely to detect the disks around any given DBZ. 
Thus, the small but non- negligible fraction of DBZ with 
observed infrared excess leads us to conclude that the 
typical disk life time is somewhat shorter than the diffu- 
sion time scales in DBZ - or that the disk lifetimes vary 
by substantial amounts (which is in fact consistent with 
the large uncertainty in our logarithmic estimate of the 
life time) . It appears unlikely that a large fraction of the 
highly polluted DBZ white dwarfs without infrared ex- 
cess are currently accreting from a gaseo us disk, as the 
life times of these gaseous disks are short (Jura 2008). 

Figure|8] also reveals a fairly well defined upper limit 
in dM/dt for DAZ white dwarfs. Based on our up- 
dated calculation of the accretion rates, there are no 
DAZ white dwarfs observed with Spitzer with accre- 
tion rates above 2 x 10 9 gs -1 . The most extremely 
metal-polluted DAZ w hite dwarf star cu rrently known, 
GALEX 1931+01170 (jVennes et afllMlfl), still only ac- 
cretes at a rate of 3 - 4 x 10 9 gs _1 (|Melis et all 120111 : 
iVennes et al.l[20TlT ). In contrast, there are several DBZ 
white dwarfs with metal accretion rates one to two or- 
ders of magnitude larger than observed in any DAZ, and 
substantially exceeding the r ates that can b e explained 
by Poynting-Robertson drag (jRafikov 2011a, tj). For in- 
stance, HE 0446-2531, HE 0449-2554, and HE 1350-162 
all have accreted several 10 23 g to nearly 10 25 g within 
the last diffusion time scale, i.e. ~ 10 6 yr (or corre- 
spondingly larger amounts of material if the accretion 
ended several diffusion time scales ago). This is compa- 
rable to the extremely metal-polluted DZ white dwarf, 
SDSS J0956+5912, w hich has at least 1.5 x 10 23 g of met- 
als in its atmosphere dKoester et al1l2011l) . Assuming the 
planetary bodies and physical mechanisms do not depend 
on the type of white dwarf, this substantial difference in 
the maximum accretion rates found for DAZ and DBZ 
white dwarfs forces us to conclude that either the calcu- 
lated metal accretion rates for DBZ white dwarfs suffer 
from systematic errors, or the extremely high metal ac- 
cretion rates found among DBZ white dwarfs simply have 
not been seen yet in DAZ white dwarfs. 

One may speculate that both DAZ and DBZ white 
dwarfs can undergo a short-lived phase of very high rate 
accretion. Whereas the long diffusion time scales of DBZ 
provide an efficient "memory" to such events, the op- 
posite is true for the DAZ, where the large amount of 
metals would be cleared out of the atmosphere within 
days to at most years - reducing the probability of wit- 
nessing such an event. A number of possible scenarios 
that would lead to such phases are conceivable. Sec- 
ondary impacts of asteroids on an existing massive dust 
disk could generate sufficient a mounts of gas to lead to 
runaway accretion (pural |2008|) . Alternatively, a colli- 
sional cascade in the initially highly eccentric disk may 
generate sufficient amounts of gas to cause runaway ac- 
cretion, effectively preventing the formation of a longer 

5 Not included in our Figure because it has not been observed 
with Spitzer 



lived disk. In either scenario, the perturbed debris disks 
are expected to evolve on relatively short timescales, and 
iGansicke et al.l (|2008l ) have shown that the structure of 
the gaseous disk in SDSS J084539. 17+225728.0 changes 
on timescales of a few years. An initial spike in the ac- 
cretion rate could also be related to the direct accretion 
of a substantial fraction of the disrupted asteroid, with 
the possible subsequent formation of a debris disk feed- 
ing the white dwarf at a lower rate. Another alternative 
comes from the fact that asteroid accretion events are, 
by definition, produced from eccentric orbits. One can 
therefore envisage that a fraction of disk are therefore 
born in an unstable configuration. The disk would, simi- 
larly to the above scenario, have many disk-disk interac- 
tions, producing a significant amount of gas. Therefore 
the disk could accreted at a very high rate, leaving a 
heavily-polluted white dwarf without a dust disk. 

5. CONCLUSIONS 

We obtained and analyzed Spitzer observations of 15 
white dwarfs with metal-polluted atmospheres, all but 
one having helium-dominated atmospheres. Of these, 
HE 0110-5630, GD61 and HE 1349-2305 are found to 
have an infrared excess consistent with closely-orbiting 
circumstellar dust. These disks are likely formed from 
the disruption of large asteroid analogs within the rem- 
nant planetary systems. A marginal excess is mea- 
sured at 15.6 /jm around NLTT 51844, but more data 
are needed to rule out contamination from extragalac- 
tic sources. 

This survey nearly doubles the number of disk detec- 
tions around DBZ white dwarfs. Using this substan- 
tially enlarged sample, we estimate a typical disk life- 
time by comparing the accreted metal masses inferred 
from DBZ stars with dust to the instantaneous accre- 
tion rates for DAZ stars with dust. Accounting for the 
large scatter by taking the logarithmic average, we find 
log[£disk(yr)] = 5.6 ± 1.1, which is compatible with the 
relatively large range of disk lifetimes estimated from dif- 
ferent theoretical models of white dwarf disks. 

The fraction of highly-polluted DBZ white dwarfs ex- 
hibiting an infrared excess is low (23 per cent) compared 
to that among DAZ white dwarfs (48 per cent). Assum- 
ing that the formation and evolution of circumstellar 
disks is similar for both types of stars, this difference sug- 
gests that the disk life times are typically shorter than 
the DBZ diffusion time scales. 

We also show that the highest time-averaged accretion 
rates are found among white dwarfs with helium-rich at- 
mospheres, many of which do not exhibit infrared excess, 
and we suggest that these stars have experienced very 
high accretion rates during short-lived phases. These 
events should occur in hydrogen-dominated white dwarfs 
as well, but their short diffusion timescales substantially 
lower the probability of detection. 
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TABLE 1 
Spitzer IRAC White Dwarf Targets 



WD 


Name 


SpT 


V 


M 


T c ff 


^cool 


r-rT /tt l 

[H/He] 


[Ca/H(e)] 


log {a. My, fat) 


log (M z ) 


Ref 








(mag) 


(M e ) 


(K) 


(Gyr) 






(gs ) 


(g) 




ni in— 565 


HE 0110—5630 


DBAZ 


15.9 





.71 


19 200 


0.12 


—4.2 


— 7.9 


8.4 


20.2 


1 2 


04.35+410 


GD61 


DBAZ 


14.9 





.72 


17300 


0.20 


—4.0 


—7.5 


8.9 


21.5 


3 


yJQ'-tO — ZOO 


I 1 1 U^l'lO — ZOOl 


FtR A 7 


lO.O 





.59 




U.OD 


1 9 


— 0. 1 


11.0 


OA Q 


1 
1 


0449—259 


T TT7* n A A n <~it~ f A 

0449—2554 


DBAZ 


16.3 





.72 


12 500 


0.52 


—4.5 


— 7.2 b 


10.0 


23.5 


1 


0802+380 


Glll-54 




15.5 





.78 


11 000 


0.85 


< —6.0 


—9.8 


7.4 


21.1 


5 


0838+375 


CBS 78 


DBZ 


17.7 





.59 


12 500 


0.36 


< -5.0 


-8.0 


9.2 


22.7 


6 


0953+594 


SDSS 


DZA 


18.4 





.58 


8 200 


1.12 


-3.2 


-7.5 


9.7 


23.6 


5 


1015+377 


CBS 127 


DZ 


17.7 





.58 


10 500 


0.59 


< -5.0 


-8.2 


9.1 


22.8 


6 


1349-230 


HE 1349-2305 


DBAZ 


16.5 





.67 


18 200 


0.14 


-4.7 


-8.0 


8.7 


22.0 


1,2 


1350-162 


HE 1350-1612 


DBAZ 


17.0 





.77 


15 000 


0.35 


-4.7 


-6.7 b 


10.4 


23.4 


i 


1352+004 


PG 


DBAZ 


15.8 





.59 


15 300 


0.20 


-4.8 


-9.3 


7.7 


20.7 


1 


1614+160 


PG 


DAZ 


15.6 





.50 


17400 


0.10 




-7.2 


7.8 




1 


2138-332 


NLTT 51844 


DZ 


14.5 





.69 


7 200 


2.23 




-8.6 


8.5 


22.5 


7 


2142-169 a 


PHL131 


DO 


15.8 


















2 


2229+139 


PG 


DBAZ 


16.0 





.83 


15 900 


0.36 


-4.5 


-9.0 


8.0 


20.9 


1 


2322+118 


LB 1188 


DZA 


16.0 





.59 


12 000 


0.41 


-5.2 


-8.7 


8.5 


22.1 


1 



Re fer ences. — (1) JKoester et al.||2005a|; (2) IVoss et al.1l2007l ; (3) IDesharnais et al.|[200Sj ; (4) IFriedrich et al.|[2000l ; (5) IDufour et~al| 
120071 ; (6) |Dupuis et al.irT993l ; ( 7) ISubasavage et al.||2007l 

Note. — Objects found to have dust disks arc listed in italics. The tenth column lists time-averaged metal in-fall rates. The eleventh 
column shows the mass of metals in the convectivc envelopes of the helium-rich stars. Both quantities arc calculated based o n the observed 
calciu m abundances and assuming this represents 1.6 per cent of the total mass of heavy elements, as in the bulk Earth [Allegre et al.| 

a Not metal-r i ch. 

k iKlein et all <|2011|) find that HS 2253+8023 has a lower calcium abundance by around 1.0 dcx than originally reported. This may 
therefore also be true of other [Fri cdrich et al. ( 2000j) stars. 
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TABLE 2 

IRAC Coordinates for HE and HS White 
Dwarfs 



Star a S 



HE 0110-5630 
HE 0446-2531 
HE 0449-2554 
HE 1349-2305 
HE 1350-1612 



01 h 12 m 21.15 s 
04 h 49 m 01.39 s 
04 h 51 m 53.72 s 
13 h 52 m 44.13 s 
13 h 53 m 34.96 s 



-56°14'27'.'8 
-25°26'36'.'l 
-25° 49' 14'.' 7 
-23°20'05'.'4 
-16°27'06'.'6 



Note. — Epoch 2009 positions as measured on 
the IRAC array from image header astrometry. 



TABLE 3 
Spitzer IRAC and IRS fluxes. 



WD 


-F3.6 M m 


-F4.5/Jm 


-FV.9Mm 






(^Jy) 


(/Jy) 


(MJy) 


(MJy) 


0110-565 


135 + 9 


124 + 8 






0435+410 




259 ± 13 


153 + 8 


98 + 10 


0446-255 




26 + 1 


15 b 




0449-259 




39 + 2 


24 + 1 




0802+386 




85 + 4 


38 + 2 




0838+375 


18 + 1 


13 + 1 






0953+594 




9 + 1 


6 b 




1015+377 


20 + 1 


14 ± 1 






1349-230 




88 + 2 


66 + 3 




1350-162 


32 + 2 


23 + 2 
23 + 1 


7 b 




1352+004 




52 + 3 


15 + 2 




1614+160 




63 + 3 


32 + 2 




2138-332 




503 ± 25 


195 ± 10 


95 + 15 


2142-169 a 




35 + 2 


12 + 4 




2229+139 




43 + 2 


13 + 2 




2322+118 


105 + 5 


69 + 3 







Note. — A 5 per cent unce rtainty has been fo lded in to 
account for calibration errors dReach et al.ll2005fl . Pipeline 
S18.7.0 was used to obtain BCD products for all targets. 
a Not metal-polluted, 
k 3<t upper limit. 
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TABLE 4 
Near-Infrared Photometry 



WD 




J 

(mag) 






H 

(mag) 






K 
(mag) 




Instrument 


0110-565 


16 


.23 


±0. 


.05 


16. 


.27 + 


.05 


16. 


24 + 


.05 


SOFI 


0435+410 


15 


.24 


±0. 


.05 


15. 


16 + 


.05 


15. 


12 ± 


.05 


LIRIS 


0446-255 


17 


.06 


±0. 


.05 


17. 


.09 ± 


.05 


17. 


15 + 


.05 


SOFI 


0449-259 


16 


.52 


±0. 


.05 


16. 


.56 + 


.05 


16. 


60 + 


.05 


SOFI 


0802+386 


15 


60 


±0. 


.05 


15. 


59 + 


.05 


15. 


65 + 


.05 


LIRIS 


0838+375 


18. 


.03 


±0. 


.05 


18. 


02 + 


.05 


17. 


.98 + 


.05 


LIRIS 


0953+594 


18 


.33 


±0. 


.06 


18. 


29 + 


.06 


18. 


23 ± 


.06 


WFCAM a 


1015+377 


17 


,24 


±0. 


.10 


17. 


31 + 


.12 


16. 


78 + 


.31 


GEMINI b 


1349-230 


16 


.91 


±0. 


.05 


16. 


94 + 


.05 


16. 


78 + 


.05 


SOFI 


1350-162 


17 


13 


±0 


05 


17. 


17 + 


.05 


17. 


.25 + 


.05 


SOFI 


1352+004 


16 


10 


±0. 


.05 


16. 


12 + 


.05 


16. 


15 + 


.05 


LIRIS 


1614+160 


15 


.93 


±0 


.05 


15. 


97 + 


.05 


16. 


.04 + 


.05 


SOFI 


2229+139 


16 


.29 


±0. 


05 


16. 


.36 + 


.05 


16. 


53 ± 


.05 


SOFI 


2322+118 


15 


.99 


±0 


.05 


15 


.98 + 


.05 


16. 


.03 + 


.05 


SOFI 



a Obtained by one of us at UKIRT using WFCAM (Casali et al. 2007, 
A&A, 467, 777). 

b Obtained by C. Melis at Lick Observatory with the GEMINI camera 
(McLean et al. 1993, SPIE, 1946, 513). 
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Fig. 1. — IRS Peak-Up image mosaics of GD61 (left) and NLTT 51844 (right). The images are orientated North up and East left at l'/8 
pixel -1 and 50" across. The background is clearly variable across the array and is therefore the primary source of error (§[27TJ. 



1.6 fim 




4.5 /xm 



7.9 fj,m 



Fig. 2. — Infrared images of GD61 taken with LIRIS at H and K s , and Spitzer IRAC at 4.5 and 7.9 /im. All images are oriented North 
up and East left, and are 30" across, with 0'.'25 pixel -1 for LIRIS and 0"6 pixel -1 for IRAC. The white dwarf (marked by the a rrow) 
is separated from a neighboring source by 1"9. The neighbor does not influence the flux measurements presented in Tablcl3l (§ 13.11 . The 
neighboring source must be in the background as its brightness never exceeds that of the white dwarf. 
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l.Q yum 



2.2 fira 



4.5 /im 



7.9 /j,m 



> 



Fig. 3. — Infrared images of HE 1349—2305 taken with SOFI at H and K, and Spitzer IRAC at 4.5 and 7.9 /im. All images are oriented 
North up and East left, and are 30" across, with 0'.'25 pixel' 1 for LIRIS and Of/6 pixel" 1 for IRAC. The white dwarf (marked by the arrow) 
is separated from a neighboring source by 2'.'8, but does not influence the flux measurements presented in Tablcl3l (§ 13. 1^ . The neighbor is 
clearly extended in the //-band image and is therefore extragalactic. 
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HE 0110-5630 




1.0 10.0 
Wavelength [yum] 



Fig. 4.— SED of HE 0110-563 0, GD61 and HE 1349-2305. The short wavelength fluxes are: HE 0110-5630: GALEX, DENIS /-band, 
and SOFI JHK, GD61: UBV ( Eggcn 1968) CMC r', and LIRIS JHK, and HE 1349-2305: GALEX, CMC r', DENIS I, SOFI JHK. 
These are shown in red with error bars. WISE data are shown in green, however, they are clearly contaminated by the neighboring source 
for HE 1349— 2305 (§3.1). The IRAC and IRS fluxes are shown in blue. The short wavelength photometry is fitted with a stellar model 
that is displayed as a black dashed line. The minimum x 2 disk model fit, shown in Figure[5] is displayed as a black dashed line, and the 
sum of the two is shown in solid black. 
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HE 0110-5630 GD61 HE 1349-2305 

90 1 — i 1 1 1 1 — i — i 1 1 1 1 — i — i 1 1 1 — 




0.3 0.7 1.1 1.5 0.3 0.7 1.1 1.5 0.3 0.7 1.1 1.5 

[1000K] 



Fig. 5.— Disk modeling for HE 0110-5630, GD61 and HE 1349-2305. The infrared excess was fitted as described in § l3~3l Briefly 
describing the format of the panels (see § 13.30 : a slice through the \ 2 volume is displayed each of the pairs of upper and lower panels for 
the minimum \ 2 fit- The upper panel shows a slice at the best fitting T ou t and the lower panel is defined by the best i. The best fitting 
solutions are: HE 0110-5630: T in = 1,000K, T out = 900 K, t = 60°, GD61: T in = 1,450K, T out = 800 K, i = 85° and HE 1349-2305: 
Ti n = 1, 700 K, Tout = 550 K, i = 85°. Regions of high \ 2 are shown as dark areas and the least \ 2 solution is marked as a red circle. Solid 
black lines show the la contours around the minimum, while the 2 and 3a contours are shown as dashed gray lines. Because only slices 
through the \ 2 are shown, this does not show the full extent of the \ 2 surface. In the case of HE 0110-5630, i vs T n , each of the regions 
enclosed by a solid black line represents a parameter space where \ 2 ls within la of the minimum. The \ 2 method does not constrain the 
disk parameters very well for HE 0110-5630. 
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Fig. 6. — SEDs of 6 of the 12 science targets consistent with photosphcric emission (§3.4). The figures follow the same general format as 
Figure|4] The SDSS optical spectrum of WD 0953+954 is shown as a gray line. The infrared excess seen in the SED of HE 0449— 2554 is 
caused by light from a background object. 
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Fig. 7. — A continuation of Figure(6]for the remaining 6 of the 12 science targets consistent with photospheric emission (§3.4). 
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Fig. 8. — Time-averaged dust accretion rates vs. cooling age of the 88 metal-polluted white dwarfs, with published abundances, observed 
with Spitzer IRAC in Cycle 1 through 7. DAZ and DBZ-type stars are plotted as open and filled circles respectively, while objects with 
infrared excess are displayed as stars rather than circles. 
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Fig. 9. — Mass of metals within the convective envelopes (or above r = 5, whichever is larger) of the 88 metal-polluted white dwarfs, 
with published abundances, observed with Spitzer IRAC in Cycle 1 through 7, plotted as a function of effective stellar temperature. DAZ 
and DBZ-type stars are plotted as open and filled circles respectively, while objects with infrared excess are displayed as stars rather than 
circles. 



